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L.  ABSTRACT 


Resxilts  of  a stnidy  of  chree  specific  items  of  incersst: 
CO  the  development  of  procedures  for  estimating  initial  radiation 
protection  factors  for  buildings  are  presented.  These  include: 

1.  The  effects  of  ne%?  cross  section  data  for 
nitrogen  and  oxygen  on  initial  radiation 
envlr'-smaents  are  presented; 

2.  Calcuia  tons  of  delayed  radiation  envixon- 
ments  ^r  large  yield  weapons  are  presimted. 

The  calculatlms  are  based  on  the  NUUSA 
code;  and 

3.  Senaltivity  analyses  of  the  effects  of 
composition,  thickness,  and  design  charac" 
terlstics  cf  wail  cons  true  tioias  are 
presented . 


I 


2.  SUMMARY 

I 

This  report  presents  the  restilts  of  a study  conducted 
for  the  Defense  Civil  Preparedness  Agency  (DCPA)  by  Science 
Applications.  Inc.  (SAI)  to  provi.de  supplementary  data  for  the 
developB»nt  of  a methodology  to  determine  initial  radiation  pro- 
tection factors  (IPF)  for  civil  defense  applications.  These  data 
involve  characterization  of  free-field  radiation  environments 
s based  on  the  beat  available  cross  section  data,  determination  of 

i the  ia^ortgnce  of  delayed  radiations  for  large  yield  nuclear 

weapons,  and  presentation  of  procedures  which  wotild  allow  the 

i 

j IPF  assessment  methodology  to  include  the  effects  of  construction 

■ material  composition  variations  and  design  characteristics. 

} Suggestions  are  presented  for  inclusion  of  the  results  of  this 

study  into  procedures  fer  estimating  IPF  for  civilian  structures. 

j 


3. 


IliTRODUCTION  AND  BACKGROUND 


Tbe  Defense  Civil  Preparedness  Agency  (DCPA)  has  for 
the  past  few  years  been  in  the  process  of  developing  procedxirea 
to  determine  the  protection  which  bnildings  will  pr'-vide  from 
initial  nuclear  radiation  (INK^  other  initial  ‘ects  asso- 
ciated with  nuclear  weapons  explosions.  These  procedures,  which 
are  similar  to  procedures  developed  for  fallout  radiation,  are 
presently  being  defined.  This  report  addresses  three  specific 
items  which  are  necessary  to  finalize  procedures  for  estimating 
initial  radiation  protection  factors  (IPF) : 

1.  The  deteimiiniation  of  the  effect  of  new  cross 
section  data  on  free  field  radiation  environ- 
ments and  their  impact  on  earlier  INR  shielding 
analysis, 

2.  The  evaluation  of  delayed  neutron  and  gamma  ray 
environments  from  large  yield  weapons,  and 

3.  The  analysis'  of  the  sensitivity  of  calculated 
dose  to  the  composition,  thickness,  and  design 
characteristics  of  walls  in  structures. 

These  three  items  are  addressed  respectively  in  Chapters  4,  5, 
and  6. 

The-,  free  field  radiation  environments  presently  being 
used  for  IPF  calculations  are  based  on  the  work  of  E.  A.  Straker 
and  M.  L.  Gritzner^^^ . These  air  transport  results  were  obtained 
in  1969  using  cross  section  data  from  ENDF/B-II.  In  the  last  few 
years,  substantial  effort  has  been  spent  by  the  Defense  Nuclear 
Agency  (DNA)  to  obtain  better  evaluations  of  the  nitrogen  and 
oxygen  cross  sections.  The  impact  of  the  newer  cross  section 
evaluations  are  reflected  in  the  analysis  described  in  Chapter  4. 
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The  development  of  the  NUIDEA  Code  sponsored  by  DNA 
and  carried  out  by  SAI  provides  a new  capability  for  estimating 
the  delayed  radiation  dose  from  large  yield  nuclear  weapons. 
Resiilts  obtained  from  this  code  are  presented  in  Chapter  5. 

For  initial  nuclear  radiation,  the  transport  of  neu- 
trons and  the  production  of  secondary  gamma  rays  is  not  only  a 
function  of  the  mass  thickness  of  shielding  material,  but  also 
a fTinction  of  the  material  composition.  Previous  calculations 
of  barrier  factors  for  neutrons  are  based  on  a fixed  concrete 
composition  and  do  not  consider  variations  in  protection  factors 
which  may  arise  from  variations  in  concrete  composition  and 
construction  methods.  Chapter  6 describes  a aeries  of  transport 
calculations  performed  for  various  thicknesses  and  several  con- 
crete compositions  of  interest.  In  addition,  perturbation  cal- 
culations are  presented  which  show  the  effect  of  specific  chemical 
constituents  of  concrete  on  the  shielding  properties.  Transport 
results  are  also  shown  for  several  wall  and  roof  constructions. 


4.  FREE  FIELD  EirrERONMENT  - PROMPT 


Previous  r'econnnendatioias  specifying  the  initial  radia- 
tions from  nuclear  weapons  were  made  in  1972  by  the  ad  hoc  Sub- 
consnittee  on  Radiation  Shielding  which  is  part  of  the  National 
Academy  of  Sciences  Advisory  Committee  on  Civil  Defense^^^ 

These  data  have  been  used  for  civil  defense  shielding  analyses 
and,  in  particular,  for  the  determination  of  radiation  protec- 
tion factors  provided  by  structures.  During  the  late  1960 's 
and  early  1970 's,  considerable  effort  was  expended  by  the 
Defense  Atomic  Support  Agency  (DASA)  and  later  the  DNA  to  develop 
methods  for  calculating  radiation  transport  in  air  and  to  improve 
the  basic  cross  section  data  reqtiired  fur  these  calculations. 

The  results  of  these  efforts  can  now  be  used  to  revise  the  free 
field  radiation  environments  recommended  in  Reference  1. 

There  are  basically  two  revisions  considered  in  this 

report.  The  first  reported  in  this  chapter  is  a revision  to  the 

prompt  neutron  and  secondary  gamma  ray  environment  due  to  updated 

(2^ 

niiclear  cross  sections  derived  from  KJA  sponsored  research-  ' . 

The  calculations  reported  in  this  chapter  were  performed  with 

(3) 

the  one-dimensional  discrete  ordinates  code,  ANISN^  ^ . The 
cross  sections  for  nitrogen  and  oxygen  used  for  the  air  trans- 
port calculations  were  obtained  from  the  DNA  few  group  library 
referenced  in  Reference  2.  Additional  problem  data  and  restilts 
of  the  calcxxlationa  are  described  below. 

4.1  TRANSPORT  PROBLEM  DATA 

4.1.1  Source  Spectra 

The  air  transport  calctilations  were  performed  using 
the  nominal  "typical”  thermonuclear  source  spectra  which  has 
been  the  convention  to  use  in  problems  of  this  type^^^ . This 
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4 


source  distributions,  grouped  in  the  DNA  few  groups  library- 

group  structxure,  is  given  in  Table  1.  Only  the  neutron  source  » 

has  been  considered.  The  omission  of  th  prompt  gamma  ray 

source  is  justifiable  because  the  pron^t  ^ ma.-’  ray  component 

of  the  prompt  dose  is  significant  only  very  near  the  source. ■ 

At  distances  of  interest  for  civil  defense  applications,  the  prompt  j 
total  dose  is  dominated  by  neutrons  and  secondary  gamma  rays. 

4.1.2  Cross  Sections  • 

The  DNA  few  group  library  is  a coupled  neutron  and 
gamma  ray  multigroup  cross  sections  library.  There  are  37 
neutron  groups  and  21  gamma  ray  groups.  The  scattering  angular 
distributions  are  approximated  by  a P2  Legendre  expansion. 

The  air  density  was  taken  to  be  0.00111  gm/cm  . At 
this  density  the  atomic  density  of  nitrogen  and  oxygen  are 
3.635x10  and  9.620x10  atoms  per  bam* cm,  respectively.  All 
other  constituents  of  air  were  conside-Jed  to  be  negligible. 

4.1.3  Response  Functions 

In  addition  to  calctilations  of  the  neutron  and  second- 
ary gamma  ray  fluxes,  tabulations  have  been  made  of  integral 
response  data  appropriate  for  estimating-  prompt  radiation  doses 
to  hunuma.  For  this  purpose  the  Snyder  Neufeld  neutron  response 
function  and  the  Henderson  tissue  gamma  ray  response  function 
have  be^  used.  This  choice  of  response  function  is  consistent 
with  that  which  has  been  used  previously.  These  response  frac- 
tions are  given  in  Tables  Z and  3. 

4.1.4  Calculational  Method 

The  neutron  and  secondary  gamma  ray  fltjx  from  a 
point  sotirce  in  air  was  calcxilated  using  the  one-dimensional 
discrete  ordinates  code,  ANISN.  An  angular  quadrattire  was 
used'.  The  transport  was  carried  to  333  go/crn^  in  air  so  that 
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Table  1.  Typical  Thcrffionuclear  Soiirce  Distribution 


Group 

Kuaber 


Upper  Energy 

Bendery  (HeV) 

Group  Fraction 
{Neutrons/Source  Neutron) 

19.6 

0. 

16.9 

0. 

14.9 

1.89(-2) 

14.2 

9. 34 (-3) 

13.8 

2. 66 (-2) 

12.8 

1.67 (-2) 

12.2 

1.69 (-2) 

11.1 

1.24 (-2) 

10.0 

7. 48 (-2) 

9. 

6.82(-3) 

8.2 

6. 78 (-3) 

7.4 

1.03(-2) 

6.4 

1.81 (-2) 

5.0 

3. 62 (-3) 

4.7 

1.24(-2) 

4.1 

2.60(-2) 

3.0 

2. 37 (-2) 

2.4 

3. 75 (-3) 

2.3 

2.56 (-2) 

1.8 

6. 44 (-2) 

1.1 

8. 85 (-2) 

5.5(-l) 

9.14(-2) 

l.S(-l) 

1.16(-2) 

1.1 (-1) 

l.ll(-l) 

5.2(-2) 

5. 40 (-2) 

2.5(-2) 

5. 68  (-3) 

2.2(-2) 

9.26(-2) 

1.0(-2) 

1.16(-1) 

3.4(-3) 

7. 38  (-2) 

1.2(-3) 

7.32{-2) 

5.8 (-4) 

2.03C-2) 

1.0(-4) 

I.90(-3) 

2.9(-5) 

0. 

l.l(-5) 

0. 

3.1(-6) 

0. 

0. 

4.7  ^-7) 
l.O(-ll) 

0. 
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Table  2.  Snyder  Naiifeld  Neutron  Response  Function. 


Group 

Huabar 

Upper  Energy 

Boundary  (MeV) 

Response  Function 
(rads /n/ cm2) 

1 

19.6 

7. 00558 (-9) 

2 

16.9 

7. 00558 (-9) 

mi 

14.9 

7. 00558 (-9) 

4 

14.2 

7. 00558 (-9) 

5 

13.8 

7. 00558 (-9) 

6 

12.8 

7. 00558 (-9) 

7 

12.2 

7.00558(-9) 

8 

11.1 

7.C0558(-9) 

9 

10.0 

7. 05279 (-9) 

9. 

7. 10289 (-9) 

U 

8.2 

7. 03619 (-9) 

12 

7.4 

6. 71089 (-9) 

13 

6.4 

6. 07429 (-9) 

14 

5.0 

5. 69619 (-9) 

15 

4.7 

5. 37649 (-9) 

16 

4.1 

4. 86219 (-9) 

17 

3.0 

4,47859(-9) 

18 

• 2.4 

4. 34239 (-9) 

19 

2.3 

4. 22839 (-9) 

20 

1.8 

3.97819(-9) 

21 

1.1 

3.34990(-9) 

22 

5.5(-l) 

1.84200 (-9) 

23 

1.6(-1) 

1.23350 (-9) 

24 

l-l(-l) 

9^.  51589  (-10) 

25 

5. 2 (-2) 

6. 92769 (-10) 

26 

2.5(-2) 

5.90470(-10) 

27 

2.2(-2) 

5.52389(-10) 

2S 

1.0(-2) 

5.57940(-10) 

29 

3.4(«3) 

6.00199(-10) 

30 

1.2<-3) 

6.16599(-10) 

31 

5.8(-4) 

6. 72759 (-10) 

32 

1.0(~4) 

5.34589 (-10) 

33 

2. 9 (-5) 

3. 88369 (-10) 

34 

l.l(-5) 

3. 43049 (-10) 

35 

3.1(-6) 

3.27479(-10) 

36 

l.l(-6) 

3.23(^0(-10) 

37 

4.1(-7) 

l.O(-ll) 

3.20529(-10) 
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Table  3.  Henderson  Tissue  Gamma  Ray  Response  Fxmction. 


Group 

Nuaber 

Opper  Energy 
Boundary  (MeV) 

Response  Function 

(rads/'':i..iac.:r.  photon/ca^) 

1 

14, 

3. 20810 (-9) 

2 

I'J. 

2.4722(-9) 

3 

8. 

2. 08470 (-9) 

4 

7. 

1.86510(-9) 

3 

6. 

1.66130(-9) 

6 

5. 

1.44310(-9) 

7 

4. 

1.19710(-9) 

8 

3- 

i.omo(-9) 

9 

2.5 

8. 70689 (-10) 

10 

2.0 

5. 64059 (-10) 

11 

1.5 

5. 64059 (-10) 

u 

1.0 

4. 10599 (-10) 

13 

0.7 

2.93009(-10) 

14 

0.45 

1.922f0(-10) 

15 

0.30 

1.10590(-10) 

16 

0.15 

5. 43209  (-U) 

17 

0.10 

3.7U30(-.U) 

18 

0.07 

3.67239(-U) 

19 

0.045 

6. 32728 (-11) 

20 

0.03 

1.41590(-10) 

21 

0.02 

0.01 

4. 40629 (-10) 
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the  spectra  will  be  free  of  botmdary  leakage  effects  to  at 
2 

least  250  gm/cm  or  a distance  of  2.25  kn. 

4.2  FSEE  FIELD  NEUTRON  AND  SECONDARY  GAMMA  RAY  SPECTRA 

Figure  la-j  shows  the  n»  .:ron  fluence  per  unit  lethargy 
per  sotiTce  neutron  for  several  ranges  from  0.2  km  to  2 km.  The 
spectra  can  be  observed  to  harden  rapidly  with  increasing  penetra" 
tion  in  air  to  about  I km.  Beyond  1 km  the  shape  of  the  spectrum 
is  changing  slowly  although  the  total  intensity  continues  to 
decrease  with  increasing  penetration. 

Figure  If  which  shows  the  fluence  per  unit  lethargy 

at  1200  ar  also  gives  a plot  of  the  energy  spectra  reconaaended 

in  Reference  1.  To  facilitate  a comparison  with  the  prompt 

environments  reported  in  Reference  1,  the  present  results  are 

also  presented  in  tht  same  format.  Figures  2a-p  show  the  free 

2 

field  neutron  fluence  multiplied  by  the  geometry  factor  4ffR 
versus  range.  Only  a few  of  the  DNA  few  group  library  groups 
correspond  directly  with  the  group  structure  utilized  in 
Reference  I. 

The  previous  liscussion  has  focused  on  spectral  dif- 
ferences which  result  by  using  more  recent  cross  section  evalua- 
tions. It  is  aL.<o  useful  to  examine  the  differences  in  the 
angular-  distributions  as  a function  of  distance  using  these 
newer  data.  A convenient  Mthod  of  illustrating  these  differ- 
ences is  to  make  direct  c(xnparison  of  the  expansion  coefficients 
of  the  angular  distributions  at  varioxis  distances  from  a point 
source.  Using,  this  procedure  provides  an  adequate  sunsnary  of 
the  isq)ortant  conclusions,  yet  presents  the  data  in  a readily 
usable  form.  Ln  Table  4,  we  present  the  ratio  of  the  harimsnic 
coefficients,  Pq  through  of  the  angular  dose  distributions 
as  a function  of  distance  from  the  source.  The  dose  response 
functions  axui  source  distributions  used  were  described  previously.  ^ 
Examination  of  these  rcstilts  indicate  that  the  changes  in  these 
distributions  becc^  more  pronotmeed  at  larger  distances , but 
are  still  quite  ssall,  l.e. , less  than  30X. 

ff 
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Figure  la.  Free  Field  Neutxon  Fluence  at  210  M. 
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Figure  lb.  Free  Field  Neutron  Fluence  at  400  H. 
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Figure  2f . Fr«e  Field  Neutron:  Group^  Fluxoe  vs  Senge. 
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Figure  2g.  Free-  Field  Neutron  Grot^  Fluxes  vs  Renge. 
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Figure  3a~j  shows  the  secondary  gamma  ray  fluence 
per  unit  lethargy  per  source  neutron  for  several  ranges  from 
0.2  to  2 km.  The  spectra  can  be  seen  to  be  relatively  insensi- 
tive to  range.  Figure  2f  compares  the  secondary  gamma  ray 
fluence  per  unit  lethargy  of  1200  m with  that  reported  in 
Reference  1.  The  spectra  are  also  plotted  as  a function  of 
range  in  Figure  4 in  the  format  of  Reference  1 to  facilitate 
a comparison  with  previous  work. 

Figures  5a-c  shows  the  Snyder  Neufeld  Tissue  dose,  the 
first  moment  of  the  dose,  and  the  second  moment  of  the  dose  as  a 
function  of  range.  Figures  6a-c  shows  the  Henderson  Tissue  gannxa 
ray  dose,  first  and  second  moments  of  the  dose,  versus  range. 
These  plots  are  normalized  to  one  so-urce  neutron. 
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5 . FREE  FIELD  ENVIRONMENT  - DELAYED 


BACKGROUND 

The  delayed  radiation  environnient  consists  of  neutrons 
and  gamma  rays  emitted  by  the  fission  products  from  a nuclear 
weapon  burst.  The  distinction  between  delayed  radiation  and 
fallout  radiation  is  that  delayed  radiation  is  that  radiation 
received  while  the  fission  products  are  contained  in  the  fireball. 
Delayed  radiation  is  defined  to  be  part  of  the  initial  radiation 
environment  which  includes  all  radiation  received  during  the 
first  minute  after  the  weapon  burst.  Even  for  yields  on  the  order 
of  10  MT,  95X  of  the  delayed  radiation  dose  is  received  within 
30  seconds  from  the  weapon  detonation. 

In  the  past,  the  pradiction  of  delayed  radiation  en- 
vironments has  been  a less  than  satisfactory  sitxaation.  The 
problem  is  a complex  one,  complicated  by  the  fact  that  the 
radiation  source  is  contained  in  a rising  and  expanding  fireball 
and  the  transport  of  radiation  is  occxirring  in  an  atmosphere  whose 
density  distribution  is  changing  due  to  the  expanding  shoe.,  front. 
For  multiple  bursts  the  problem  is  further  complicated  by  inter- 
acting shock  fronts  and  the  possibility  of  interacting  fireballs 
and  is  not  well  xmderstood.  Although,  in  principle,  the  problem 
cotild  be  solved  using  radiation  hydrodynamics  codes,  the  physical 
size  of  the  problem  precludes  this  as  a practical  meai\s  to  a 
solution.  Historically,  empirical  models  based  on  test  data  ha.va 
been  developed  and  applied  to  the  prediction  of  the  delayed  en- 
vironments. Unfortunately,  sufficient  data  are  not  available 
« 

to  construct  an  empirical  model  for  the  range  of  yields,  burst 
heights,  and  weapon  types  of  interest,  ifore  recently  the  problem 
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has  been  approached  by  phenomenology  modeling  as  exemplified 
by  Che  computer  code,  NUIDEA^^\ 

THE  NUIDEA  CODE 

The  development  of  the  NUIDEA  Code  was  sponsored  by 
the  Defense  Nuclear  Agency  (DNA) . The  documentation  and  code 
shoxxld  be  available  to  users  in  the  defense  conanunity  in  the 
near  future. 

The  NUIDEA  Code  was  developed  as  a systems* like  code 
for  the  investigation  of  nuclaar  weapon  environments  from  single 
or  multiple  bursts.  The  code  includes  phenomenology  models  of 
nuclear  radiation,  blast,  and  thermal  radiation.  The  code  in- 
corporates  portions  of  the  Air  Transport  of  Radiation  (ATR)  '■ 

Code  and  the  Low  Altitude  Multiple  Burst  (LAMB) Code. 

5.1  CALCULATIONS 

A matrix  of  calcvilations  were  performed  using  the  NUIDEA 
Code  for  weapon  yields  frcm  1 KT  to  10  MT  and  for  three  burst  heights 
in  meters,  1 meter,  60  m,  and  225  m (where  W is  the  yield 

in  kilo  tons) . The  scaled  biirst  heights  were  selected  as  chose  burst 
heights  which  precliwie  fallout  and  optimize  blast  effects  (>10  psi) , 
respectively.  Several,  components  of  the  delayed  radiation  environ- 
ment were  tabulated  including  tissue  dose  for  delayed  gamma  rays, 
delayed  neutrons,  and  secondary  gamma  rays  from  delayed  neutrons. 

The  protapt.  radiation  dose,  blast  overpressure,  and  thermal  exposure 
are  also  provided  by  the  NUIDEA  Code  and  have  been  tabulated  as 
well.  For  each  waapon  yield  considered,  the  fission  yield  was  as- 
sumed to  be  80%  of  the  total  yield.  The  ground  elevations  for  these 

cases  was  380  meters. 

* 

5.2  RESULTS 

5.2.1  Radiation  Envirodtaents 

Figures  7,  8,  and  9^  show  the  radiation  eavlronssnts  at 
UOO  Rt  ground  versus  yield  for  the  three  burst  heights, 


ground  range 

HOB  • 60w^/3 

'^OTAL  DOSg 

total 

prompt  dose 

oeuY  neutron 


OELA?  2NO 


YIEU)  (Jcn 


FlgtJre  S. 


S«>S«  f“  S0*W  K 


respectively.  The  figiires  show  all  the  components  of  the  ini- 
tial radiation  dose  including  both  prompt  and  delayed  coi^onents. 
It  is  observed  for  the  groimd  burst  case  that  for  yields  greater 
than  a few  htmdred  kilotons  the  delayed  gamma  dose  dominates  the 
total  dose.  However,  for  the  60  case  the  delayed  gamma  dose 

exceeds  the  total  prompt  dose  only  for  yields  greater  than  2 MT. 
For  the  225  case,  the  prompt  dose  is  always  the  predominate 

component  of  the  total  dose. 

Figures  10,  11,  and  12  show  a similar  set  of  curves 
for  a ground  range  of  2500  meters.  The  san»  observations  apply 
as  in  the  1500  meter  ground  range  cases. 

5.2.2  Blast  Overpressure 

Figure  13  ahows  the  maximum  overpresstnre  from  a 1 MT 
burst  versus  ground  range  for  the  three  burst  heights.  The 
ground  range  for  a given  overpressure  can  be  deteioained  for 
other  yields  based  on  the  standard  cube  root  scaling  laws 

^^ere  GR^  is  the  ground  range  at  yield 

GRq  is  the  ground  range  at  yield  Wq 

5.2.3  Thermal  Exposxire 

Figures  14,  15,  and  16  show  the  total  thermal  expostjre 
versus  gromd.  range  and  weapon  yield  for  the  three  respective 
burst  heights. 
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6.  TRANSPORT  THSOTIGH  STRUCTURES 


This  chapter  presents  results  of  calculations  of  the 
transport  of  neutrons  and  gansaa  rays  through  structural  materials. 
The  transport  of  neutrons  include  the  prodtiction  of  secondary 
gamma  rays  and  their  transport.  Results  for  three  types  of  cal- 
culations are  presented 

1.  One-dimensional  slab  transport  calculations, 

2.  Sensitivity  calculations  for  transport  through 
concrete,  and 

3.  Ring  sourca  effects. 

The  structures  considered  in  the  Iculations  include 
» Concrete  slabs, 

• Wood  frame  walls, 

• Brick  veneer  walls, 

• ShJaigle  roofs,  and 

• Bxiilt  up,  asphalt  roofs. 

Forward  and  .adjoint  one-dijnensional  transport  calculations  were 
performed  using  the  ANISN^^^  discrete  ordinates  code.  The  for- 
ward calculations  used  a source  distribution  determined  directly 
from  tiw  free  field  calculations  reported  in  Ch.apter  4.  All  cal- 
culations used  cross  sections  from  the  DNA  few  group  library. 


6.1 


TRANSPORT  THROUGH  CONCRETE 


6.1.1  Concrete  Compositions 

Concrete  does  not  uniquely  specify  an  elemental  com- 
position in  the  stnse  required  for  radiation  transport  analysis, 
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Variations  in  moisture  content,  regional  variations  in  the 
conjtituents  of  concrete,  and  the  basic  type  of  concrete  all 
will  effect  the  shielding  properties.  Our  approach,  therefore, 
was  to  select  several  basic  types  of  concrete  for  analysis  and 
to  account  for  variations  within  a particular  type  by  perturba- 
tion methods. 

The  elemental  composition  for  several  important  types 
of  concrete  are  listed  in  Table  5. 


6.1.2 


Soxarce  and  Response  Functions 


The  sotirce  used  for  the  forward  one-dimensional  slab 


transport  calculations  was  a "shell  source"  taken  directly  from 
the  air  transport  calculations  reported  in  Chapter  4.  The 
ANISN  Code  provides  the  capability  of  coupling  calculations  in 
this  manner.  The  assumptions  made  in  coupling  the  calculations 
in  this  manner  is  (1)  the  flux  incident  on  the  slab  is  not  per- 
turbed by  the  presence  of  the  slab  and  (2)  the  radius  of  the 
coupling  sxorface  is  sufficiently  large  and  lateral  transport 
sufficiently  limited  that  the  switch  from  spherical  to  slab 
geometry  is  appropriate. 

The  response  fxjnction  which  was  also  used  as  the 
soizrce  distribution  for  the  adjoint  calculations  was  the 
Snyder  Neufeld  neutron  response  and  the  Henderson  Tissue  gamma 
ray  response  fisaction  (refer  back  to  Chapter  4,  Section  4.1.3 
and  Tables  2 and  3) . 


6.1.3 


Forward  Transport  Results 


The  total  dose  transmission  as  a function  of  areal 
mass  is  plotted  in  Figxire  17  for  the  seven  types,  of  concrete 
listed  in  Table  4.  Figure  18  shows  the  ratio  of  the  neutnron 
dose  to  the  gaimna  ray  dose  aa  a function  of  areal  mass.  Both 
the  transmission  and  the  n/y  ratio  can  be  seen  to  vary  sub- 
stantially with  concrete  type  e^'evi  at  the  same  areal  mass. 
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Table  5 . Concre 


StilaJdtB)  coM(K>«l 


DOSE  THANS«^IS$)ON 


MKirraoN  oose/tjAWMA  hav  ooss  hatjo 


6.L.4 


Perturbation  Calculations 


In  order  to  investigate  the  sensitivity  of  the  trans-* 
mitted  dose  through  concrete  to  the  concrete  composition,  a 
series  of  perturbation  calculations  were  performed  tising  the 
SAIDOT  Code^^^ . The  sensitivity  calculation  is  based  on  the 
perturbation  relationship 


AR 

X 


(»'^,AL») 

R 


where 


R is  the  transmitted  dose  for  a reference  problem, 

♦ is  the  flux  for  the  reference  problem, 

is  the  adjoint  flux  for  the  reference  problem, 

AL  is  the  "perturbation"  to  the  transport  operator,  and 
AR  is  the  "perturbation"  to  the  transmitted  dose. 


In  the  present  calctilations  the  perturbation  to  the  transport 
operator,  AL,  has  been  considered  to  be  a change  in  the  concrete 
composition.  By  using  the  perturbation  relationship  and  having 
calculated  the  forward  and  adjoint  fluxes  for  a reference  problem, 
the  transmitted  dose  for  any  concrete  composition  can  be  estimated 
as  long  as  the  difference  in  the  concrete  composition  is  not  too 
different  from  the  composition  of  the  reference  problem.  If  we 
define  the  sensitivity  function  for  a particxilar  element  to  be 


where 

then 


where 


is  the  microscopic  cross  section  for  element  i. 

# elements 

R'  - RfAR  - R[l+^(N£-Nj^)S^] 

i^l 

K*.  is  the  atonic  density  of  eleaent  t in  the  reference 
concrete  and 

is  the  atomic  dr-nsity  of  elenent  t in  the  "perturbed" 

case. 
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Table  6 gives  the  sensitivity  ftmction  as  a func- 
tion of  concrete  thickness  using  type  TSF  concrete  as  a reference 
case. 

The  range  of  applicability  for  the  pertnirbation  can  be 
tested  by  using  the  perturbation  results  to  predict  the  trans- 
mission factors  for  the  concrete  coo^ositions  listed  in  Table  5 
and  comparing  vd.th  the  transport  results.  Table  7 shows  the 
fractional  error  of  the  permrbation  theory  prediction  when  com- 
pared with  the  transport  (ANISN)  results.  Using  TSF  concrete  as 
a base  case,  perturbation  theory  can  be  seen  to  accurately  pre- 
dict the  transmitted  dose  for  ordinary  concrete  Type  03  and  Type 
04,  as  well  as  for  Serpentine  concrete.  However,  pert\arbation 
theory  does  not  give  an  accurate  prediction  for  the  Magnetite 
concrete,  and  Limonite  with  steel  punchings.  Perturbation  theory 
fails  for  these  cases  because  the  Iron  loading  for  these  con- 
cretes is  so  high.  However,  it  appears  that  pertxirbations  from 
TSF  concrete  can  be  quite  accurately  predicted  for  concrete  con- 
taining less  than  about  5 weight  percent  iron.  In  order  to 
treat  the  more  heavily  iron  loaded  concretes  with  perturbation 
theory  would  require  a reference  fluxes  for  a more  heavily 
loaded  concrete.  The  restilts  based  on  TSF  concrete  shovild  be 
acctirate  for  most  concretes  cotsnonly  encountered  in  the  con- 
struction industry. 

6.2  TBANSPOKT  THROUGH  OTHER  STRUCTURAL  ELEMENTS 

The  transmission  factors  for  other  structural  elements 
commmly  found  in  the  building  industry  have  also  been  calculated. 
These  structural  elements  include  roof  and  wall  constructions 
('.ooBonly  fotmd  in  residential  hoiBes.  These  include  both  ^od 
And  brick  exterior  walls,  the  shingle  roof,  and  the  built  up 
asphalt  and  gravel  roof.  Figtnre  19  shows  the  materials  and  con- 
figuration of  these  stmictural  elements.  In  order  for  these 
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configuracions  to  be  amenable  to  one-dioKnsional  tranaport 
analyses,  the  wood  studs  were  approximated  by  a hooKSgenotis 
region  of  1/8  the  density  of  wood  (the  ratio  2:16).  The  re- 
stilts  of  the  one- dimens tonal  transport  calculations  for  these 
configurations  are  summarized  in  Table  8. 

6.3  RING  SOURCE  EFFECTS 

The  calculations  reported  in  Sections  6.1  and  6.2 
for  dose  transmission  through  concrete  slabs  and  other  build- 
ing structural  elements  have  been  for  a point  source  above  an 
infinite  slab.  Since  adjoint  calculations  were  performed  for 
the  concrete  slabs  it  is  particularly  convenient  to  fold  these 
adjoint  fluxes  with  the  incident  fltixes  from  other  sotirce  con- 
figurations. The  ring  source  is  of  partictxlar  Interest  since 
it  has  been  the  recommended  source  configuration  for  civil 
defense  shielding  analysis.  The  ring  source  results  presented 
here  are  based  on  a rotation  of  the  point  sotirce  fluxes  incident 
on  the  slab  followed  by  a convolution  with  the  slab  response 
function. 

6.3.1  Source  Rotation  and  Lengendre  S?.pansion 

The  problem  is  to  calculate  the  flux  on  the  axis  of 
a ring  source  in  air,  given  the  flux  from  a point  source  in  air. 
This  can  be  determined  by  a straightforward  rotation  of  coordin- 
ates. Let  [$1  represent  a vector  whose  elements  are  the  flux 
mon^nts  for  the  point  source. 
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Table  8.  RestdCs  of  Walls  and  Roofs  Calculation. 


Configuration 

Uall: 

Wood  Exterior 
Brick  Exterior 

Roof: 

Shingle 
Built  Op 


Doae, 

rad 

source  neutron 

i 

Neutron 

Stay 

Total 

Dose 

Transalsalon 

3.6«(-22) 

1.46(-22) 

5. 14 (-22) 

0.37 

2.46(-22) 

1.01 (-22) 

3. 48 (-22) 

0.25 

2.86<-22) 

1.48(-22) 

4. 34 (-22) 

0.32 

2,28(-22) 

1.50<-22) 

3. 78 (-22) 

0.28 

where 


4'  = J $(u)Pj^(u)du  , 

*•  L 

u » cosine  of  the  angle  measured  from  the 
line  to  the  source 

are  the  Legendre  polynominals . 

Also,  let  represent  the  corresponding  Legendre  expansion 
of  the  flux  for  the  ring  source,  then 


where 

and 


where 


- M[«] 

4 - « 

JPj^Ccos  e)dj 
^ -1  0 

e and  are  illustrated  below. 


RING  SOURCE 


a can  be  determined  from  8q,  9,  and  4 as  follows. 

a. 

u = cos  0Q  cos  6 + sin  bQ  sin  0 cos  0 . 

6.3.2  Ring  Source  Restilts 

The  source  rotation  and  folding  with  the  adjoint  flux 
was  perforts»d  for  three  thicknesses  of  TSF  concrete  and  for  ring 
source  delimation  angles.  9q,  from  0 to  90®  (0®  corresponds  to 
the  point  sotnrce) . The  ratio  of  the  transmitted  dose  from  the 
ring  soiirce  to  the  transmitted  dose  from  the  point  sovrrce  are 
plotted  in  Figure  20  as  a function  of  the  cosine  of  the  decline- 
tiOTi  angle.  In  general,  the  transmitted  dose  decreases  with 
increasing  declination  angle  due  to  the  effectively  Increased 
a-wtrage  path  length  through  the  slab.  It  is  interesting  to  note, 
however,  that  for  the  thinner  slabs,  build  up  effects  cause  a 
small  increase  in  the  transmitted  dose  over  the  first  few  degrees 
of  declination.  These  results  used  the  infinite  air  angular 
fluxes  at  1200  m. 


^^■54  gm/csn^ 


M.88  g^cnr^ 


COS^ 


Figure  20. 


Ring  So\irce  Effects  for 
Sc-irce  at  1200  H in  Air 


a TherssKjnucXear 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  relative  importance  of  delayed  radiation  for  civil 
defense  applications  for  high  yield  weapons  has  been  demonstrated. 
These  data  are  presented  in  Chapter  5.  The  angular  dependence  of 
the  delayed  radiation  have  not  been  studied  in  detail,  however, 
preliminary  indications  are  that  the  angular  distributions  tend 
to  be  more  forward  peaked  then  the  prompt  radiation  and,  there- 
fore, could  have  some  impact  on  the  wall  and  roof  barrier  factors 
for  INR.  These  effects  should  be  evaluated. 

The  use  of  newer  cross  section  data  to  determine  the 
free  field  environments  from  prompt  radiation  indicates  some 
differences  relative  to  the  ENDF/B-II  data.  These  differences 
should  be  indicated  in  the  final  methodology  which  is  to  be  used 
for  INR.  applications. 

It  appears  that  the  selective  use  of  perttirbation  tech- 
niques is  more  than  sufficient  to  determine  that  variation  (either 
an  increase  or  decrease)  in  the  initial  protection  factor  for 
changes  in  material  compositions  based  on  elemental  differences . 

It  is  recommended  that  provisions  be  made  in  the  INR  methodology 
for  inclusion  of  procedures  to  estimate  the  effects  of  both  ele- 
mental composition  and  construction  techniques  on  predicted  values 
of  IPF.  These  techniques  could  quite  easily  be  developed  by  using 
the  data  in  Cliapter  6,  along-  with  .some  supplementary  calculations. 
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